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This paper presents the modeling, design, and simulation of an adaptive 
neuro fuzzy inference strategy (ANFIS) for controlling the speed of the 
Double Star induction Machine (DSIM), the machine is fed by three phase 
direct matrix converter which makes directly AC-AC power conversion is 
modeled using Direct Space Vector Modulation technique(DSVM) for direct 
matrix converter. Double star Induction motor is characterized by highly 
non-linear, complex and time-varying dynamics and inaccessibility of some 
of the states and outputs for measurements. Hence it can be considered as a 
challenging engineering problem in the industrial sector. Various advanced 
control techniques has been devised by various researchers across the world. 
Some of them are based on the neuro-fuzzy techniques. The main advantage 


Neuro-fuzzy of designing the ANFIS coordination scheme is to control the speed of the 

Space vector modulation DSIM to increase the dynamic performance, to provide good stabilization. 
To show the effectiveness of our scheme, the proposed method was 
simulated on an electrical system composed of a 4.5 kW six-phase induction 
machine and its power inverter. Digital simulation results demonstrate that 
the deigned ANFIS speed controller realize a good dynamic of the DSIM, a 
perfect speed tracking with no overshoot, give better performance and 
high robustness. 
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1. INTRODUCTION 

The induction motor takes a large place in industry and we can find it in electrical ships. However, 
when high reliability is required, the control using the field-oriented method is often necessary and on-line 
diagnosis or off-line diagnoses have to be considered. The use of a double star induction motor or six-phase 
induction motor is slowly increasing and we find it in the high power process. Its main advantage lies in most 
reliability in case of a normal operating system [1, 2]. Double star machine supplied with non sinusoidal 
waveforms causes perturbations in the torque (harmonics in torque). The induction motor drive fed by a 
matrix converter is superior to the conventional PWM-VS inverter because of the lack of bulky DC-link 
capacitors with limited life time, the bi-directional power flow capability, the sinusoidal input/output 
currents, and adjustable input power factor. Various modulation methods for matrix converters have been 
investigated. Space Vector Modulation (SVM) is one of the preferred modulation methods for three -phase 
matrix converter [3]. The Space Vector Modulation (SVM) techniques are the extension of the theory of flux 
in multi-phase rotating machines to the field of static power converters. SVM strategy enables the control of 





Journal homepage: http://beei.org/index.php/EEI 


Bulletin of Electr Eng and Inf ISSN: 2302-9285 0 819 





input power factor independent from output power factor, achieves maximum voltage transfer ratio without 
the need of third harmonic injection and reduce harmonics [4]. 

The DSIM it is desirable to control the flux and torque separately in order to have the same 
performances as those of DC motors. One way of doing this is by using the field oriented control. 
This method assures the decoupling of flux and torque. The vector-controlled DSIM with a conventional PI 
speed controller is used extensively in industry, because has easily implemented. Alongside this success, 
the problem of tuning PlI-controllers has remained an active research area. Furthermore, with changes in 
system dynamics and variations in operating points PI-Controllers should be returned on a regular basis. 
Today there are many methods for designing intelligent controllers, such as predictive controller, 
fuzzy control, neural networks and expert systems. Various combinations of these controllers give a number 
of design possibilities. Artificial Neural Networks (ANNs) and Fuzzy Logic (FL) have been increasingly in 
use in many engineering fields since their introduction as mathematical aids. A combination of neural 
networks and fuzzy logic offers the possibility of solving tuning problems and design difficulties of fuzzy 
logic. The resulting network will be more transparent and can be easily recognized in the form of fuzzy logic 
control rules or semantics [5]. This new approach combines the well established advantages of both the 
methods and avoids the drawbacks of both. 


2. INDIRECT FIELD ORIENTED CONTROL OF DOUBLE INDUCTION MACHINE 

The model of the machine that is established is described by a set of differential equations linking 
statoric and rotoric magnitudes. The dynamic behavior of this machine supplied by two current inverter is 
described by a model with rotor flux as state variable [6]. 


del, 1 Ms" ; 
08 = 2 oh + agegh +E Cita + iha) () 
deh 1 Mg" ‘ 
cae SS ar Pq + WgMq + ae (ig + bq) (2) 


@g: being the slip pulsation. 
The electromagnetique torque is described by the next equation 


T= 2p: (ME) [oh(ita + a) — 8 Ca + a] ~ 


P1: number of the pairs vectorial control by orientation of the flux necessitates a choice of judicious 
referential. The choice of a particular referential allows transforming the expression of the electromagnetic 
torque so that the induction machine can be compared to D.C machine at least by torque expression. In our 
study, we have chosen the case of an indirect control by orientation of the rotor flux, whose advantages 
resides in the simplification and the decoupling of equations of the system. In the model of the double-star 
asynchronous machine represented by (1) to (3), we choose a referential linked to the rotating field in order 
that the direct axis D coincides with the desired direction of the rotor flux [7]. This choice of referential 
imposes ~p"(c) = gG(c) = 0 and @G(c) = 0. Thus (1) to (3) can be written as: 


det i . 
TS + 04 = ME (ita + ida) (4) 
_ ME Gigtizg 
Og = To (5) 
3 Msr ‘ . 3 oh 
T. = 5Pa (“) palitg + iba) = 5 P1 ar Og (6) 


The mechanical equation is given by: 


ja" + £0" =T, —T, (7) 
With: T.: torque of load 

j: moment of inertia 

f: damping constant of the machine 
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The model that is established constitutes the main idea of the control of the double star 
asynchronous machine. We notice that components if ,andi5,. If the rotor flux, while the torque depends only 
of quadrature components of the stator currents i?,andi3,. If the rotor flux g’ is maintained constant. 
For the double star asynchronous machine supplied by current inverter. The stator currents (i7g, (3q, i{q, tq) 
and the slip pulsation are considered as variables of control. 


3. MATRIX CONVERTER MODELING AND CONTROL 
3.1. Modeling for matrix converter 

A matrix converter is a variable amplitude and frequency power supply that converts the three phase 
line voltage directly. It is very simple in structure and has powerful controllability. The real development of 
the matrix converter starts with the work of Venturini and Alesina who proposed a mathematical analysis and 
introduced the low frequency modulation matrix concept to describe the low frequency behavior of the 
matrix converter [8]. In this, the output voltages are obtained by multiplication of the modulation matrix or 
transfer matrix with the input voltages. The basic diagram of a matrix converter can be represented 
by Figure 1. The existence function provides a mathematical expression for describing switching patterns. 
The existence function for a single switch assumes a value of unity when the switch is closed and zero when 
the switch is open. 





Figure 1. Basic structure of matrix converter 


For the matrix converter shown in Figure 2, the existence function for each of the switches is expressed by 
the following equations: 


1, switch § kj closed 
ki = (8) 


0, switch S,; open 


Where k={A, B, C} is input phase and j={a, b, c} is output phase. The above constraint can be expressed in 
the following form: 


Sa +83 +Sq =1 (9) 


with the above restrictions a 3x3 matrix converter has 27 possible switching states. The mathematical 
expression that represents the operation of a three phase ac to ac Matrix Converter can be expressed 
as follows: 


v(t) Sg) Spa) Sol) vy(t) 
vp) ]=|S4,0 SppQ) Sq@)|*| vg (10) 
v.(t) Sac) Spe) So Vc(t) 
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in| [Sua Spa Sad] [ig 
ing(t)|=|Sap@) Spot) S(t] *]i,0 (1) 
ic(t) S4-(t) S p(t) So, (t) i.(t) 


where va, vb and ve and iA, iB and iC are the output voltages and input currents respectively. To determine 
the behavior of the MC at output frequencies well below the switching frequency, a modulation duty cycle 
can be defined for each switch. The modulation duty cycle MKj for the switch SKj in Figure 2 


is defined as in (12). 


t.: 
M,=— (12) 


where tkj is the one time for the switch Skj between input phase k={A, B, C} and j={a, b, c} and Ts is the 
period of the PWM switching signal or sampling period. In terms of the modulation duty sycle, 
(13)-(14) can be rewritten as given below. 


v(t) Mygt) Mz) Ma va(t) 
vpt)|=|M 4,0 Mgy(t) Mey) |*| ve(t) (13) 
v.(t) M,.(t) Mz.) Me@)} lve) 


ig] [Maalt) Malt) Ma] [ia 
ip(t)}|=|My,) Mz) Ma) */ i, 
ic) | |My) Mz.) Me.) i,(t) 
My,+My+MgG =1 (14) 
j={a, b, c} 
The high-frequency synthesis technique introduced by Venturini and Alesina in [9, 10] allows the 
use of low frequency continuous functions, referred to as the modulation matrix m(t), to calculate the 


existence functions for each switch of the matrix converter. Thus, the aim when using the Alesina and 
Venturini modulation method is to find a modulation matrix which satisfies the following set of equations. 


v,(t) = m(t).v;(t) (15) 
i) =m)’ i, (16) 
where the input voltages vi(t) are given by the following set of functions 
cos(@;t) 
2. 
v,(t) = V,,| cos(@;t - = (17) 
cos(@,t + 285 
3 
and the desired output voltages vo(t) are 
cos(@,t) 
v, (t) =V,| cos(@,t — =) (18) 


cos@,t + =) 
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output currents i, (¢) can be expressed as: 


cos(@,t + @,) 
i,(@) =T1, c0s(oyt- + py) (19) 


2a 
cos(@,t + 3 +Q,) 


Where 0 is the phase angle of the linear load. Finally, the desired input current has an arbitrary phase qi. 
This angle can be set to 0 to obtain unity input power factor of the matrix converter. 


cos(@,t + Q;) 
it) =I; cos(eat ~ == +9.) (20) 


2a 
COS(@,t + a +@Q;) 
The elements of matrix m(t) that satisfy (15) and (16) are given by 


m,(t) = sali + 4004, —@;)+ SH - |} + 
(21) 


2 aft geo o, ~0)+200-1-p)} 
31 3 


V, 
+———— |, @, =1- qa 


where, a, =— 
tan(@, ) i 


1 I tan(Q, ) 
2 


With the following restrictions a, 20,a,20,0<q< 


Nl 


3.2. Direct space vector modulation technique for matrix converter 

The space vector modulation (SVM) is a control technique that has been widely used in adjustable 
speed drives and more generally in power converter control. The idea of Direct Space Vector Modulation 
(DSVM) for the matrix converter, unified representation of the current and the voltage space vectors is 
presented [11]. In an output of reference frame, the instantaneous voltage space vector of the output 
is given by (23). In an input of reference frame, the instantaneous current space vector is given by (24). 


B= =(an + Vay + FVeq) = vi(theli© (22) 
Te = = (vag + F vac + Gvcq) = Vo (tele (23) 
T= =(ig tT ip + Fig) = i(eIO (24) 
ig = = (ig t+ 7 ig + Pic) = ig (teHPo (25) 


where vi(t), vo(t) , 1i(t) and io(t) are the time dependent magnitudes of the space vectors while ai(t), ao(t), 
Bi(t) and Bo(t) are the corresponding time dependent phase angles. For a three phase MC there are 27 valid 
switch combinations giving thus 27 voltage vectors as shown in Table 1. The switching combinations can be 
classified into three groups which are, synchronously rotating vectors, stationary vectors 
and zero vectors. 

The DSVM algorithm selects the switching states and calculates the duty cycle for each 
switching state. During one switching period Ts, the switching states are adjacent to the desired output 


Bulletin of Electr Eng and Inf, Vol. 8, No. 3, September 2019 : 818 — 828 


Bulletin of Electr Eng and Inf ISSN: 2302-9285 0 823 





voltage (input current) vector, should be selected, and the zero switching states are applied to complete the 
switching period to provide maximum output to input voltage transfer ratio. Figure 2 shows input and output 
hexagons, and Table 2 presents the MC states that have to be chosen in terms of reference vectors location. 





Table 1. Possible SCS and vectors used in the MC 





states Switches on |Z, | (@, liol G, 
ABB +1 Sas Sav Spe 2/3 Vap 0 2/V3 ig 2/6 
BAA -1 Spa Sab Sac —2/3 Vag 0 -2/V3 ig —n/6 
BCC 42 Spa Sch Sco 2/3 Vac 0 2/V3 ia n/2 
CBB -2 Sca Sap Spe 2/3 Vgc 0 -2/V3 ig 0/2 
CAA 43 Sca Sav Sac 2/3 Vea 0 2/V3 ig 71/6 
ACC -3 Saa Scp Sco 2/3 Veg 0 -2/V3 ig  7n/6 
BAB +4 Spa Sav Spe 2/3 Vag 20/3 2/V3 i, 2/6 
ABA -4 Saa Sps Sac  —2/3 Vag 20/3 -2/V3 i, 2/6 
CBC +5 Sca Spb Sco 2/3 Vac 20/3 2/V3 ip n/2 
BCB -5 Spa Scp Spe —2/3 Vgc 21/3 —2/v3 ip m/2 
ACA +6 Saa Sco Sac 2/3 Ve, 20/3 2/Vv3 ly 7n/6 
CAC -6 Sca Sav Sco 2/3 Veq 20/3 -2/V3 ip 7n/6 
BBA +7 Spo Spo Sac 2/3 Vag 40/3 2/V3 i,  —n/6 
AAB -7 Si Se Sas —2/3 Vag 41/3 -2/V3 i, —1/6 
CCB +8 Sca Sco Spe 2/3 Vac 41/3 2/Vv3 ic 1/2 
BBC -8 Spa Spb Sce —2/3 Vac 41/3 —2/V3 ic 1/2 
AAC +9 Saa Sav Sce 2/3 Vea 41/3 2/v3 ig 7n/6 
CCA -9 Sca Seo Sac 2/3 Vea 4r0/3 -2/V3 i,  7n/6 
AAA O, Saa Sav Sac 0 0 
BBB 0 Spa Spy Spe 0 0 





Table 2. Active vector configurations 
Output voltage reference vector sector 








lor4 2 or 5 3 or 6 
Input current lor4 9 7 3 1 6 64 9 7 3 1 6 4 
reference vector 2or5 8 9 2 3 5 6 8 9 2 3 5 6 
location 3 or 6 - 8 1 2 4 5 7 8 1 2 4 5 
I TI @wm ImWvI. I I IV I II Il IV 





The duty cycle of each selected state id introduced to denote the on-time ratio of that particular 
switching state during one switching period Ts. 51, . 62. 63 and . 64 are defined as the duty cycles of the four 
switching states I, II, III and IV respectively. Kv and Ki denote the sector in which the output voltage and 
input current space vectors are located in the af reference frame. The duty cycles are calculated based on the 
phase of output voltage and input current vector references such as in the following [12]: 


2m 05(06-*)cos(a-2 


6, = —1hot Kina as <o (26) 
sage ™ 
ede a ea 0 
64 = —iket hint am cos(9o+5)eos( +4) (29) 


v3 cos(A@) 


where m is the modulation index: Ag is the displacement angle between the measured input voltage vector vi 
and the input current reference vector ii*, kv and ki are the voltage and current sector, and: 


5 =D — (ky -— DF (30) 


0, = 0; - (ki; -1)> (1) 
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m = ¥3/2 


The duty cycles must satisfy (32) in order that the short circuit of any two input phases or open circuit of the 
load terminals does not occur. 


15,| + ]52] + 163] + [64] < 1 (32) 
The duty cycle for the zero vector is calculated using (33). 


yo = 1— {16,] + 162| + [53] + 184} (33) 


+4+526 





5 
+3+6i9 








(a) (b) 


Figure 2. Available Vectors for SVM, (a) Voltage vectors, (b) Current vectors 


4. DESIGN OF ANFIS CONTROLLER 

The acronym ANFIS derives its name from adaptive neuro-fuzzy inference system. Using a given 
input/output data set, the toolbox function anfis constructs a fuzzy inference system (FIS) whose membership 
function parameters are tuned (adjusted) using either a backpropagation algorithm alone, or in combination 
with a least squares type of method. This allows your fuzzy systems to learn from the data they are modeling. 

ANFIS’s network organizes two parts like fuzzy systems. The first part is the antecedent part and 
the second part is the conclusion part, which are connected to each other by rules in network form. If ANFIS 
in network structure is shown, that is demonstrated in five layers. It can be described as a multi-layered 
neural network as shown in Figure 3. Where, the first layer executes a fuzzification process, the second layer 
executes the fuzzy AND of the antecedent part of the fuzzy rules, the third layer normalizes the membership 
functions (MFs), the fourth layer executes the consequent part of the fuzzy rules, and finally the last layer 
computes the output of fuzzy system by summing up the outputs of layer fourth. [13, 14]. 

In order to model complex nonlinear systems, the ANFIS model carries out input space partitioning 
that splits the input space into many local regions from which simple local models (linear functions or even 
adjustable coefficients) are employed. The ANFIS uses fuzzy MFs for splitting each input dimension; the 
input space is covered by MFs with overlapping that means several local regions can be activated 
simultaneously by a single input. As simple local models are adopted in ANFIS model, the ANFIS 
approximation ability will depend on the resolution of the input space partitioning, which is determined 
by the number of MFs in ANFIS and the number of layers. Usually MFs are used as bell-shaped with 
maximum equal to | and minimum equal to 0, The feed forward equations of ANFIS are as follow [15, 16]: 


[y)(x)=——+— (34) 


b+ (—7y" 


i 


Mai(2) = exp {- yh (35) 


a; 
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Where {ai, bi, ci} are the parameters of MFs which are affected in shape of MFs. In this study the ANFIS 
controller generates change in the reference voltage Vier, based on speed error e and derivate in the speed 
error de defined as: 

€ = Wref — W (36) 

de = [d(@;er — w)]/dt (37) 


where ref and @ are the reference and the actual speeds, respectively. Figure 3 shows the ANFIS 
model structure. 


input inputmt rule outputmt output 


Logical Operations 
and 


@ « 


ms not 





Figure 3. The ANFIS model structure 


5. SIMULATION RESULTS 

The SIMULINK model for indirect FOC of the 4.5 Kw cage rotor DSIM associated with ANFIS 
controller is shown in Figure 4. The machine is fed by a matrix converter with Space vector modulation 
Technique. The parameters of the induction motor are summarized in Appendix. The first test concerns 
a no-load starting of the motor with a reference speed @,.¢ = 288 rad/sec. and a nominal load disturbance 


torque (14N.m) is suddenly applied between Isec and 2sec, followed by a consign inversion (-288rad/sec) 
at 2.5sec. At 4.5s, a -14Nm load disturbance is applied during a period of 2 s. this test has for object the study 
of controller behaviors in pursuit and in regulation. The test results obtained are shown in Figure 5. 

The speed of the motor reaches @,,¢ at 0.2 s with almost no overshoot. It then begins to oscillate 


inside a 0.4% error strip around ®,.¢ . The ANFIS controller rejects the load disturbance very quickly with no 


overshoot and with a negligible steady state error. In order to test the robustness of the used method we have 
studied the effect of the parameters uncertainties on the performances of the speed control. To show the 
effect of the parameters uncertainties, we have simulated the system with different values of the parameter 
considered and compared to nominal value (real value). The Figure 6 and Figure 7 show respectively the 
behavior of the DSIM when Rr is 10% increased of its nominal value and J is increased 10% of its 
nominal value. An increase of the moment of inertia gives best performances, but it presents a slow dynamic 
response. The figures show that the proposed controller gave satisfactory performances thus judges that the 
controller is robust. 
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Figure 4. Simulink diagram for DSIM control systems 
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Figure 5. Simulated results of ANFIS controller for DSIM 
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Figure 6. Simulated results of ANFIS controller for DSIM with variation 
of the rotor inertia (+10%J) 
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Figure 7. Simulated results of ANFIS controller for DSIM with variation 
of the rotor resistance at t=2s 


6. CNCLUSION 

This paper presents a study of an application of speed control with neuro fuzzy inference strategy 
(ANFIS) controller for a double stator induction machine based on the direct FOC, The machine is fed by 
a matrix converter, Simulation results on control robustness with speed variation parameters variation and 
the torque resistant are given. The results show that the controller could compensate for this kind of 
disturbances. The plant is also tested for the tracking property using different types of reference signals. 
Satisfactory performance was observed for most reference tracks and the results demonstrated the 
effectiveness of the proposed structure and the proposed control scheme it is believed will constitute a major 
step in the evolution of intelligent control 
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APPENDIX 


Table 3. DSIM parameters 








N nominal values Value IS-Unit 
Power 4.5 kW 
Frequency 50 Hz 
Voltage ( A TY) 220/380 Vv 
Current ( A TY) 5.6 A 
Speed 2751 rpm 
Constant Value IS-Unit 
Poles of pair 1 
Rs = Rs 3.72 Q 
Rr 2.12 Q 
Ls,= Ls) 0.022 H 
Lr 0.006 H 
Lin 0.3672 H 
J 0.0625 Kgm? 
Ky 0.001 Nm(rad/s)! 
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